Spitzer Planet Limits around the Pulsating White Dwarf GD66 
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ABSTRACT 

We present infrared observations in search of a planet around the white dwarf, GD66. Time- 
series photometry of GD66 shows a variation in the arrival time of stellar pulsations consistent 
with the presence of a planet with mass >2.4Mj. Any such planet is too close to the star to 
be resolved, but the planet's light can be directly detected as an excess flux at 4.5 /im. We 
observed GD66 with the two shorter wavelength channels of IRAC on Spitzer but did not find 
strong evidence of a companion, placing an upper limit of 5-7 Mj on the mass of the companion, 
assuming an age of 1.2-1.7Gyr. 

Subject headings: infrared: stars — planetary systems — stars: oscillations — white dwarfs 



1. Introduction 

The advantage of searching for sub-stellar com- 
panions around white dwarf stars (WDs) has long 
been recognized. The improved contrast between 
a WD and a companion was first exploited by 
IProbstI (|l983l ) and the first L dwarf around a WD 
(or indeed any star) was d i scover ed around GD165 



bv iBecklin fc Zuckerman (|1988f ). In recent years, 
improved instrumentation has enabled searches 
for lower mass objects, to the point where di- 
rect detection searches for planetar y mass objects 



are poss ible from the gr ound (e.g. iBurleigh et al 
[2002: Debes et all l2005l Hogan et al. 2008, in 



press). Searches utilizing the Spitzer Space Tele- 
scope (jWerner et al. 20041) can ob tain mass lim- 
its as low as ^ 10 Mj (IFarihi et al.ll2008a : .MuUallvl 



20071: iDebes et al.ll2007l) 



WDs are interesting targets for more reasons 
than just their luminosity. By looking at WDs 
we can sample a range of progenitor masses from 
about 1-8 M0, including higher mass stars not 
amenable to radial velocity surveys. Currently, 
the best main-sequence targets for direct detec- 
tion of planetary companions are either transiting 
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syste ms ( Charbonneau et al. 2005t Deming et al 
I2OO6 , e.g.) or very young stars ( Marois et al.ll2008t 



Ka las et al.l 120081 ). Detection of a planet around 



a WD would open up a older, more sedate re- 
gion of parameter space where the predictions of 
planet models could be tested against observation 
for ages ranging from hundreds of millions to bil- 
lions of years. WD planets are key to testing our 
models of planetary evolution. 

WDs also lend insi ght into the u l timate fate 



of planetary systems. iLivio fc Sokerj (|1984| ) first 



estimated the conditions necessary for a gas gi- 
ant planet to survive the red g iant p hase of its 
parent star. ISackmann et al.l ( 1993f ) concluded 



that the Earth would survive the red giant phase 
of the Sun, but more r e cent calculations by 



Schroder fc Connon SmithI (|2008l ) concluded the 



opposite. These calculations depend on stellar 
mass loss rates that are weakly constrained by 
observation due to the rapid evolution of evolved 
stars. Determining the distribution of planets 
around WDs will provide a key test for these 
models. The winds from evolved stars pollute 
the interstellar medium with the material that 
creates the next generation of terrestrial planets, 
so surveys for planets around dead stars will in- 
form our understanding of planets around stars 
not yet born. 

Early empirical success comes from the detec- 
tion of a planet around a sub-dwarf star using pul- 
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sation timings (|Silvotti et al 



2007). planets in an 



eclipsing sub-dwarf system ()Lee et al. 20081 ) and 
the detect ion of silicate rich debris disks aroun d 
WDs (e.g. iReach et al.ll2005bHFarihi et al.ll2007l ). 
Analysis of the composition of material raining 
down onto the sur f ace of a WD from such a disk by 
IZuckerman et al. ( 2007 ) concluded that the pat- 
tern of elemental abundances is consistent with a 
species of asteroid. Asteroids are markers of plan- 
etary systems, and so the frequency of debris disk 
WDs may eventually provide an independent mea- 
sure of the fraction of stars that bear (or bore) 
planetary systems. 

iMuUallv et alJ I l|2008 l. hereafter MG8) published 
the results of a survey that sought to detect evi- 
dence of planets as variations in the arrival time of 
pulsations from variable WDs (in a manner anal- 
ogous to sub-dwarf and pulsar planets). For one 
object, GD66, they discovered a variation consis- 
tent with a 2.1 Mj planet in a 4.5 yr orbit at a 
distance of about 50 pc. Unfortunately, the span 
of their data was not long enough to cover an en- 
tire orbit and their detection remains provisional. 
In this paper we report on follow up observations 
with Spitzer in an attempt to confirm this object 
by direct detection in the infrared. Although the 
candidate system is too distant to resolve into indi- 
vidual components, the low luminosity of the star 
means that the flux from the planet at 4.5 fim ri- 
vals that of the WD. For wavelengths less than 
10 /imatmosphere models of planets are brightest 
at 4-5 /im, in a gap betw een absorb ption bands of 
methane and water (e.g iBurrows et al. 2003). The 
flux from this bump increases with mass, and de- 
creases with age, and provides a method of both 
detecting a planet and determining its mass. The 
goal of these observations is to measure an excess 
flux at 4.5 /xm over the bare photosphere of GD66 
due to the light from the planet. We achieve this 
by measuring the flux ratio (i.e. the color) between 
chann els 2 and 1 on IRA C (4.5 and 3.6 /xm respec- 
tively, iFazio et"aL 2004 ) and comparing this ratio 
to a sample of other similar stars. This approach 
is capable of detecting companions independent of 
distance to the star, or the assumptions of atmo- 
sphere models. 

These observations push the IRAC instrument 
in a new direction. While studies of transit- 
ing planets with IRAC have been able to detect 
changes in flux in a single band of order a few 



parts in 10~^, our goal was to measure an absolute 
flux excess relative to other wavelengths. With 
this approach, we established a sensitivity limit of 
~0.5% with the two IRAC bands available during 
the warm mission. 

2. Observations 

We observed GD66 with all four IRAC chan- 
nels on 17th October 2007 (AORs 22855424 and 
22856448) as part of program 40255. We recorded 
sixteen 100 second exposures simultaneously in 
channels 1 and 3, using a spiral dither pattern. 
We repeated the procedure for channels 2 and 4, 
recording 4 sets of sixteen images for a total ex- 
posure time of 6,400 seconds. The raw frames 
were calibrated using version S16.10.0 of the IRAC 
pipeline to create bed frames. 

We also observed 3 WDs with similar tem- 
peratures to calibrate our observations. L19- 
2 (AORs 22856192 and 22857216), ZZ Ceti 
(AORs 22855936 and 22856960) and G238-53 
(AO Rs 22855680 a n d 228 56704) are all DAVs, 
and iBergeron et al. ( 2004f ) reports temperatures 
within 120 K of GD66 based on optical spec- 
troscopy. This agreement is better than the 
quoted uncertainty in the measurements, and so 
we can consider our sample to have a uniform 
temperature. As an independent, albeit cruder 
approach, we also insisted on bro adly similar pul 



sation properties for each star, as lMukadam et al 



(l2006l) showed that pulsation properties are cor- 
related with temperature. We observed each of 
these stars in a similar pattern as GD66, scaling 
our total exposure time with magnitude. 

We extracted photometry from these frames us- 
ing the astrolib package in IDL. We flrst convert 
the recorded flux into units of electrons. Pixels 
with bits set in any of the pmask, imask or rmask 
images were marked as bad and their values set 
to the median value of pixels in a box 21 pixel on 
a side centered on the bad pixel in question. We 
ignored bit 13 of the imask, indicating crosstalk, 
which was set in a large number of pixels that did 
not seem to be affect ed by any noticeable p roblem. 
Graham et all (|l990l ) fc lReach et all (|2008h see ev- 
idence of pulsation at infrared wavelengths in G29- 
38, a WD of similar tempe rature but also pos - 
sessing a close-in debris disk (jReach et al.ll2005bf ) . 
However, we find no evidence of pulsation or any 
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systematic drift in the measured flux in any chan- 
nel for any star in our sample. 

We extract the flux using a sky annulus of 10- 
20 pixels and series of apertures ranging from 2 to 
7 pixels radius. We choose the aperture that gives 
the greatest signal to noise (2 pixels), adding in 
quadrature a contribution to the error term due to 
the read noise, as listed in the header of the bed 
frame. We then apply corrections to the observed 
flux as suggested by the Data Handbook: array lo- 
cation dependence, aperture correction, color cor- 
rection and a pixel phase correction for Channel 
1. 

3. Results 

In Figure [1] we show spectral energy distribu- 
tions (SEDs) for GD66, L19-2 and ZZ Ceti. We 
show a table of observed fluxes for all observed 
stars in Table [1] and the flux ratio between chan- 
nels 2 and 1, r, for each star in Figure [2] and Ta- 
ble H 

GD66 and G238-53 are noticeably redder 
(higher value of r) than the other two stars. How- 
ever, careful analysis of the images shows that this 
flux excess is most likely an artefact for G238-53. 
The diffraction spike from a nearby bright star 
passes very close to the position of G238-53 on 
the chip. Unlike the other three stars, the obser- 
vations for the two channels for this star were not 
take sequentially, but separated by ~ 1 week. The 
slightly different viewing angle between the two 
visits means that the diffraction spike is closer to 
the star in the epoch of channels 2 and 4 (which 
are observed simultaneously). Although some of 
these pixels are masked, we consider the measured 
fluxes to be untrustworthy. The flux ratio is sig- 
nificantly higher (0.675 versus 0.657) for a 3 pixel 
aperture, a further indication that flux is con- 
taminated. We therefore ignore G238-53 in our 
analysis. 

3.1. Error Analysis 

We now determine whether this excess in r for 
GD66 is statistically significant. The error budget 
of an IRAC observation is divided to photomet- 
ric, systematic and calibration components. The 
IRAC data handbook gives a value of 3% for the 
calibration uncertainty while Farihi et al. ( 2008bh 
concludes, based in part on a literature survey 



that 5% is more appropriate. However, because 
we are measuring a flux ratio, we need not con- 
cern ourselves with the absolute calibration error 
in converting from electrons to mJy necessary to 
compare IRAC fluxes to other photometric sys- 
tems. We do, however, include a conservative un- 
certainty term of 5% in Table [T] to facilitate com- 
parison with other stars. 

The photometric uncertainty in channel 1 can 
be computed from the rms scatter of individual 
bed frames. Adding the fractional uncertainties 
in quadrature gives the fractional uncertainty in 
the flux ratio. Our observations were designed to 
achieve <1% photometric accuracy, and this goal 
was comfortably exceeded, as shown in Tabled] 

The systematic uncertainty is measure of the 
variability in the performance of the instrument, 
and the spread in values we would obtain by per- 
forming the same experiment multiple times. The 
Observers' Manual quotes a stability for IRAC 
of <1%. To more accurately measure this value, 
we examined data on BD -1-60 1753, a calibration 
star that has been repeatedly obser ved by IRAC 



over the lifetime of the instrument ( Reach et al 



2005al ). We reduced 81 AORs on this star span- 
ning a period of approximately 3.6 years with the 
same algorithm as we used for our science targets. 
For each AOR we measured the average of 5 point- 
ings as before, and divided the flux at 4.5 /im by 
that at 3.6 /xm. We show a histogram of this ra- 
tio in Figure [21 The typical photometric error bar 
on these measurements is 5r=0.0023, and the rms 
deviation is 0.0042. As our error terms add in 
quadrature, the systematic repeatability limit is 
0.0035, or 0.54%. 

We add our independent estimate of the sys- 
tematic uncertainty in quadrature to the photo- 
metric uncertainty to get a total (1 cr) uncertainty 
in the ratio, Sr, which is given in the third column 
of Table [2] and shown in Figure [21 The excess flux, 
E is the difference between r for GD66, 0.6482(62) 
and the mean value of r for the other two stars, 
0.6365(31), or £:==0.0117(69). The excess flux is 
only 1.7 times greater than the uncertainty and 
therefore not statistically significant. 

4. Limits on Planet Mass 

While we did not detect an excess, this does not 
rule out the evidence from timing measurements, 
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and we can use our observations to constrain the 
mass of any companion. We can state with 3a con- 
fidence that the flux ratio is less than 0.6556. To 
compare this number to planet models and place 
an upper bound on the mass of the companion we 
first need the age of the system. 

Based on the measured surface temperature 



and g ravity (11,989 K and log q— 8.05 lBergerpn et al 
2004 ). and comparing to models from Iwood 



(|1992I ). we determine a mass of 0.64(03) and 
a cooling age of 500 Myr. The age of the pro- 
genitor main-sequence star is more difficult to 
estimate. Estimates of the initial-final mass re- 
lationship (IFMR) for WDs are based on com- 
paring the masses of WDs in open clusters to 
the main-sequence turn-off age of the cluster it- 
self. Until recently, this relationship was only 
empirically constrained for young, and hence high 
mass progenitors. However, newer observations 
have begun to extend the IFMR down to the rela- 
tively lower ma ss of GD66. Using the linear fits to 
the IF MR f romlDobbie et al.l (l2006f ) , iKalirai et al 
fioO?) and 'Catalan et al.' (l2008l) we estimate an 
initial mass of 2.20(46), 2.26(46) and 2.64(57) Mg 
respectively, which agree within th e uncertain- 
ties. A theoretical relationship from iMeng et~al 



(|2007n gives an intermediate value of 2.49(41) Mp 



while the semi-empirical method of iSalaris et al 
( 20081) su ggests 2.31(22) M©. Calculations from 
Pols et al . (1998) indicate that the lifetime of a 
2.5 M0 star is 830 Myr. The lifetime of a main se- 
quence star scales roughly as M~^-^ in this mass 
regime, so the calculated main-sequence lifetime 
of GD66 is 700 Myr to 1.2 Gyr. Including the 
WD cooling age, the total age of the system is 
1.2-1.7 Gyr. 

We take atmosphere models of a 12,000 K DA 
courtesy of D. Koester (jFinlev et al ][l993), and 
combine t hem with models of p lanetary mass ob- 
jects from lBurrows et aL ( 2003f ). We perform syn- 
thetic photometry on the combined model to de- 
termine the flux ratio between channels 1 and 2 
for planets of different masses. We repeat the pro- 
cedure for models with ages from 0.3 to 5 Gyr and 
show the results in Figure ID Given our limit on 
the flux excess, and the uncertainty in the age, 
the 3(T upper bound on the mass of the planet is 
approximately 5-6 Mj . 



5. Additional Time Series Observations 

We also present additional ground-based opti- 
cal tim e-series observations of GD 66 taken with 
Argos (|Nather fc MukadamI |2004| ) on the 2.1m 
telescope at McDonald Observatory. A full de- 
scription of these observations is given in M08, but 
we repeat a brief summary here. We performed 
time series differential photometry on GD66 with 
a continual sequence of exposures over a period of 
4-8 hours per night. We divided the lightcurve of 
the target star by the sum of one of more reference 
stars to remove the effects of seeing and trans- 
parency variations and we corrected our times for 
the motio n of the earth a round the sun using the 
method of Stumpfi ( 1980l ). We combined all obser- 
vations over the span of a week into a single data 
set, and measured the difference between the ob- 
served time of arrival of the 302s pulsation mode 
and that calculated based on the assumption of a 
constant period (0-C). We show our updated 0-C 
diagram in Figure [5l 

M08 predicted that the orbit would turn over in 
late 2007 and the slope of the 0-C diagram would 
be negative in 2008. This has clearly not hap- 
pened and the orbital period is longer than pre- 
viously expected. Our current best fit parameters 
are a period of 5.69(30) yr, an orbital separation 
of 2.75(11) AU and mp.smi = 2.36(17) Mj, where 
TOp is the mass of the planet and i the inclination 
of the orbit to the line of sight. With previous ex- 
perience in mind we are cautious in predicting the 
future behavior of the orbit, but we note that our 
current circular best fit orbit appears to be on the 
verge of turning over. 

6. Discussion 

We did not detect a planet around GD66, but 
our observations showcase the advantages of tar- 
geting WDs in direct detection planet searches. 
Stars at the end point of their evolution are many 
magnitudes fainter than when on the main se- 
quence, which more than compensates for the de- 
clining fiux from the cooling planet. This im- 
proved constrast, together with the typically high 
proper motions of WDs make them ideal candi- 
dates for surveys with ground based high contrast 
imager s being built for lar ge telescopes (e.g. Hi- 
CIAO, iTamura et a"Lll2006f ). 
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6.1. Excesses at Long Wavelengths 
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We observe slight excesses in channel 4 for L19- 
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Table 1 
Observed fluxes 



Star 


J 


Integration 


Ch 1 


Ch 2 


Ch3 


Ch4 




(mag) 


(s) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


GD66 


15.7 


6400 


0.1349(68) 


0.0875(44) 


0.0593(50) 


0.0277(35) 


L19-2 


13.6 


768 


0.941(47) 


0.600(30) 


0.383(22) 


0.219(12) 


ZZCeti 


14.4 


1920 


0.470(24) 


0.299(15) 


0.189(12) 


0.1105(67) 


G238-53 


15.7 


6400 


0.1385(69) 


0.0911(46) 


0.0573(39) 


0.0405(26) 



Note. — Integration times are for channels 2 and 4, and are four times shorter for 
channels 1 and 3. The flux from G238— 53 is probably overestimated due to the presence 
of a diffraction spike from a nearby star. The uncertainty in the two least significant 
digits is given in parentheses. The quoted uncertainties in these measurements include 
a 5% absolute calibration uncertainty added in quadrature. 



Table 2 

Observed fluxes ratios and uncertainties 



Star 


Ratio 


Photometric 


Systematic 


Total 




(Ch2/Chl) 


Uncertainty 


Uncertainty 


Uncertainty 


GD66 


0.6482 


0.0051 


0.0035 


0.0062 


ZZCeti 


0.6364 


0.0026 


0.0034 


0.0043 


L19-2 


0.6376 


0.0028 


0.0034 


0.0044 


G238-53 


0.6574 


0.0030 


0.0036 


0.0046 
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Fig. 1. — Fluxes of three WDs in all 4 IRAC bands. The sohd line connecting the points is a 12,000 K DA 
white dwarf atmosphere model. The model is normalised to agree with the flux at 3.6/L(mand is shown to 
guide to the eye. 
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Fig. 2. — Flux ratio between channel 2 and channel 1 for three WDs. The dash line is the weighted average 
flux of the two comparison WDs. 
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Fig. 3. — Histogram of ratio of measured flux between channel 2 and channel 1 for the IRAC calibration 
star, BD +60 1753 over a period of 3 years. While the photometric scatter in any one AOR on this object 
predicts an uncertainty in this ratio of about 0.002, the measured rms scatter is about 0.004. This systematic 
difference is primarily due to long term drifts in the instrumental response. See text for details. 
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Fig. 4. — Mass limits on a planetary companion as a function of the assumed age. The solid lines show the 
expected flux ratio for WD+planet systems of different ages. The shaded region indicates the best estimate 
of the age of GD66. The downward arrow indicates the 3a upper bound on the flux ratio. 
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Fig. 5. — 0-C diagram for GD66. Filled circles are the observed arrival times of the pulsations. Values 
greater than zero indicate arrival times that were later than expected based on the assumption of a constant 
pulsation period and no companions. The solid line is the best fit sinusoid to the data. 
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